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Abstract 

We use the Maple system to check the investigations of S. S. Gupta regarding the Smaran- 
dache consecutive and the reversed Smarandache sequences of triangular numbers [Smaran- 
dache Notions Journal, Vol. 14, 2004, pp. 366-368]. Furthermore, we extend previous 
investigations to the mirror and symmetric Smarandache sequences of triangular numbers. 
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The nth triangular number tn, n G N, is defined by tn = X^ILi ^ = n{n + l)/2. These numbers 
were first studied by the Pythagoreans. 

The first k terms of the triangular sequence {tn}^=i are easily obtained in Maple: 

> t:=n->n*(n+l)/2: 

> first := k -> seq(t (n) ,n=l . . k) : 

For example: 

> first (20) ; 

1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 66, 78, 91, 105, 120, 136, 153, 171, 190, 210 

In this short note we are interested in studying Smarandache sequences of triangular numbers 
with the help of the Maple system. 

To define the Smarandache sequences, it is convenient to introduce first the concatenation 
operation. Given two positive integer numbers n and m, the concatenation operation cone is 
defined in Maple by the following function: 

> cone := (n,m) -> n*10~length(m)+m: 
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For example, 

> coiic(12,345) ; 

12345 

Given a positive integer sequence {un}^=i, we define the corresponding Smarandache Con- 
secutive Sequence {scsn}'^=i recursively: 

SCSI = Ui , 

SCS„ = COnc{sCSn-l,Un) ■ 

In Maple we define: 

> scs_n := (u,n) -> if n=l then u(l) else conc(scs_n(u,n-l) ,u(n)) fi: 

> scs := (u,n) -> seq(scs_n(u, i) , i=l . .n) : 

The standard Smarandache consecutive sequence, introduced by the the Romanian mathe- 
matician Florentin Smarandache, is obtained when one chooses Un = n,'in E N. The first 10 
terms are: 

> scs(n->n, 10) ; 

1, 12, 123, 1234, 12345, 123456, 1234567, 12345678, 123456789, 12345678910 

Another example of a Smarandache consecutive sequence is the Smarandache consecutive 
sequence of triangular numbers. With our Maple definitions, the first 10 terms of such sequence 
are obtained with the following command: 

> scs(t,10) ; 

1, 13, 136, 13610, 1361015, 136101521, 13610152128, 
1361015212836, 136101521283645, 13610152128364555 

Sometimes, it is preferred to display Smarandache sequences in "triangular form" . 

> show := L -> map(i->print(i) ,L) : 

> show([scs(t,10)]) : 

1 

13 
136 
13610 

1361015 
136101521 
13610152128 
1361015212836 
136101521283645 
13610152128364555 
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The Reversed Smarandache Sequence (rss) associated with a given sequence {un}'^=i, is 
defined recursively by 

rssi = ui , 

rssn = conc{un, rssn-i) ■ 
In Maple we propose the fohowing definitions: 

> rss_n := (u,n) -> if n=l then u(l) else conc(u(n) ,rss_n(u,n-l)) fi: 

> rss := (u,n) -> seq(rss_n(u, i) , i=l . .n) : 

The first terms of the reversed Smarandache sequence of triangular numbers are now easily 
obtained: 

> rss(t,10) ; 

1, 31, 631, 10631, 1510631, 211510631, 28211510631, 
3628211510631, 453628211510631, 55453628211510631 

We define the Smarandache Mirror Sequence (sms) as follows: 

smsi = ui , 

smsn = conc{conc{un,smSn-i),Un) 

> sms_n := (u,n) -> if n=l then 

> u(l) 

> else 

> cone (cone (u(n) , snis_n(u,n-l) ) ,u(n)) 

> fi: 

> sms := (u,n) -> seq(sms_n(u, i) , i=l . .n) : 

The first 10 terms of the Smarandache mirror sequence introduced by Smarandache are: 

> sms (n->n, 10) ; 

1, 212, 32123, 4321234, 543212345, 65432123456, 7654321234567, 
876543212345678, 98765432123456789, 109876543212345678910 

We are interested in the Smarandache mirror sequence of triangular numbers. The first 10 
terms are: 

> sms (t, 10) ; 

1, 313, 63136, 106313610, 1510631361015, 21151063136101521, 
282115106313610152128,3628211510631361015212836, 
45362821151063136101521283645,554536282115106313610152128364555 

Finally, we define the Smarandache Symmetric Sequence (sss). For that we introduce the 
function "But Last Digit" (bid): 
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> bid := n -> iquo(n,10): 

> bld(123); 

12 

If the integer number is a one-digit number, then function bid returns zero: 

> bld(3); 



This is important: with our cone function, the concatenation of zero with a positive integer 
n gives n. 

> concCbldd) ,3) ; 

3 

The Smarandache Symmetric Sequence (sss) is now easily defined, appeahng to the Smaran- 
dache consecutive, and reversed Smarandache sequences: 

SSS2n-l = COnc{bld{sCS2n-l),rSS2n-l) , 
SSS2n = COnc{sCS2n, rSS2n) , 

n G N. In Maple, we give the following definitions: 

> sss_n := (u,n) -> if type(n,odd) then 

> conc(bld(scs_n(u, (n+l)/2)) ,rss_n(u, (n+l)/2)) 

> else 

> conc(scs_n(u,n/2) ,rss_n(u,n/2)) 

> fi: 

> sss := (u,n) -> seq(sss_n(u, i) , i=l . .n) : 

The first terms of Smarandache's symmetric sequence are 

> sss(n->n, 10) ; 

1, 11, 121, 1221, 12321, 123321, 1234321, 12344321, 123454321, 1234554321 
while the first 10 terms of the Smarandache symmetric sequence of triangular numbers are 

> sss(t,10) ; 

1, 11, 131, 1331, 13631, 136631, 136110631, 1361010631, 1361011510631, 13610151510631 

One interesting question is to find prime numbers in the above defined Smarandache sequences 
of triangular numbers. We will restrict our search to the first 1000 terms of each sequence. 
All computations were done with Maple 9 running on a 2.00Ghz Pentium 4 with 256Mb RAM. 
We begin by collecting four lists with the first 1000 terms of the consecutive, reversed, mirror, 
and symmetric Smarandache sequences of triangular numbers: 

> st:=time(): LscslOOO: = [scs(t, 1000)] : printfC'/.a seconds" .round (time ()-st)) ; 
20 seconds 

> st:=time(): LrsslOOO : = [rss (t , 1000)] : printfC'/.a seconds" .round (time ()-st)) ; 
75 seconds 
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> st:=time(): LsmslOOO: = [sms(t,1000)] : printf ("%a seconds" ,roviiid(time()-st) ) ; 



212 seconds 



> st:=time(): LssslOOO: = [sss(t,1000)] : printf("7.a seconds" .round (time ()-st)) ; 
26 seconds 

We note that scsiooo and rssiooo are positive integer numbers with 5354 digits; 

> length(LscslOOO[1000]) , length (Lrss 1000 [1000] ) ; 

5354, 5354 

while smsiooo and sssiooo have, respectively, 10707 and 4708 digits. 

> length(LsmslOOO[1000]) , length (Lsss 1000 [1000] ) ; 

10707,4708 

There exist two primes (13 and 136101521) among the first 1000 terms of the Smarandache 
consecutive sequence of triangular numbers; 

> st : = time () : 

> select (isprime jLscslOOO) ; 

> printf ("Za minutes" ,round( (time () -st) /60) ) ; 

[13, 136101521] 

9 minutes 

six primes among the first 1000 terms of the reversed Smarandache sequence of triangular 
numbers; 

> st : = time : 

> select (isprime .LrsslOOO) ; 

> printf ("Za minutes" ,round((time()-st)/60)) ; 

[31, 631, 10631, 55453628211510631, 786655453628211510631, 10591786655453628211510631] 

31 minutes 

only one prime (313) among the first 600 terms of the Smarandache mirror sequence of 
triangular numbers;^ 

> length(Lsmsl000[600]) ; # sms_{600} is a number with 5907 digits 

5907 

> st : = time () : 

> select (isprime,Lsmsl000[l. .600]) ; 

> printf ("Za minutes" ,round( (timeO -st) /60) ) ; 

[313] 

^Our computer runs low in memory when one tries the first 1000 terms of the Smarandache mirror sequence 
of triangular numbers. For this reason, we have considered here only the first 600 terms of the sequence. 
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3 minutes 



and five primes among the first 1000 terms of the Smarandache symmetric sequence of trian- 
gular numbers (the fifth prime is an integer with 336 digits). 

> St : = time : 

> select(isprime,LssslOOO) ; 

> printf("%a minutes" ,round((time()-st)/60)) ; 

[11, 131, 136110631, 1361015212836455566789110512012010591786655453628211510631, 
1361015212836455566789110512013615317119021023125327630032535137840643546549652856159 
5630666703741780820861903946990103510811128117612251275132613781431148515401596165316 
5315961540148514311378132612751225117611281081103599094690386182078074170366663059556 
152849646543540637835132530027625323121019017115313612010591786655453628211510631] 

19 minutes 

> lengthen [5]); 

336 

How many primes are there in the above defined Smarandache sequences of triangular num- 
bers? This seems to be an open question. 

Another interesting question is to find triangular numbers in the Smarandache sequences of 
triangular numbers. We begin by defining in Maple the boolean function istriangular. 

> istriangular := n -> evalb(nops(select(i->evalb(whattype(i)=integer) , 

> [solve (t(k)=n)])) > 0): 

There exist two triangular numbers (1 and 136) among the first 1000 terms of the Smarandache 
consecutive sequence of triangular numbers; 

> St := timeO : 

> select (istriangular jLscslOOO) ; 

> printfCZa seconds" ,round(time()-st) ) ; 

[1, 136] 

6 seconds 

while the other Smarandache sequences of triangular numbers only show, among the first 
1000 terms, the trivial triangular number 1: 

> St : = time : 

> select (istriangular ,LrsslOOO) ; 

> printf("°/oa seconds" ,round(time()-st) ) ; 

[1] 

6 seconds 
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> St := timeO : 

> select (istriangular, LsmslOOO) ; 

> printf ("°/oa seconds" .round (time ()-st) ) ; 

[1] 

10 seconds 

> st : = time : 

> select (istriangular, LssslOOO) ; 

> printf ("7oa seconds" ,round(time()-st) ) ; 

[1] 

6 seconds 

Does exist more triangular numbers in the Smarandache sequences of triangular numbers? 
This is, to the best of our knowledge, an open question needing further investigations. Since 
checking if a number is triangular is much faster than to check if a number is prime, we invite 
the readers to continue our search of triangular numbers for besides the 1000th term of the 
Smarandache sequences of triangular numbers. We look forward to readers discoveries. 
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